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ABSTRACT 

How important is an independent diameter measurement for the determination of stellar parameters of solar-type 
stars? When coupled with seismic observables, how well can we determine the stellar mass? If we can determine 
the radius of the star to between 1% and 4%, how does this affect the theoretical uncertainties? Interferometry 
can provide an independent radius determination and it has been suggested that we should expect at least a 4% 
precision on such a measurement for nearby solar-type stars. This study aims to provide both qualitative and 
quantitive answers to these questions for a star such as our Sun, where seismic information will be available. 
We show that the importance of an independent radius measurement depends on the combination of observables 
available and the size of the measurement errors. It is important for determining all stellar parameters and in 
particular the mass, where a good radius measurement can even allow us to determine the mass with a precision 
better than 2%. Our results also show that measuring the small frequency separation Sh' significantly improves the 
determination of the evolutionary stage r and the mixing-length parameter a. 

Subject headings: stars: oscillations — stars: fundamental parameters (mass, age, initial hydrogen abundance, 
initial metal abundance, mixing-length parameter) — stars: interiors — methods: numerical 
— techniques: interferometric 



1. INTRODUCTION 

Asteroseismology is the interpretation of a star's oscillation 
frequency spectrum to characterize its internal structure. By 
probing its interior, we are testing our knowledge of fundamen- 
tal physics. A star with oscillations such as our Sun (solar- 
type star) presents a frequency spectrum with a range of excited 
modes. Because most of the proper modes are excited, mode 
identification for each frequency i s considerably easier than 
for other types of pulsating stars (Kield sen & Bedding 120041: 
[Bedding et al. 2004; Kjeldsen et al. 2005). 

Modelling these stars can be difficult as there are often mul- 
tiple sets of parameters that fit the observations. Also, if we use 
the frequencies to determine the mass of the star, then any solu- 
tion we find may be incorrect as it depends on the^f mass. If we 
could determine the mass of the star through some other means, 
then modelling the frequencies would be easier, and they could 
also then be used solely to characterize the internal structure of 
the star 

Unfortunately the precise determination of mass is normally 
possible only for detached components of similar mass in spec- 
troscopic binary systems. For single stars we use observ- 
ables such as magnitudes and colors (apart from oscillation 
frequencies) to estimate the mass. However, now that there 



are various possibilites of measuring a radius for si ngle stars 
through interferometry (Paresce et all 2003: Glindem ann et al.l 
l2003l:lKervella et alj 120031: lAufdenberg'et a l. 2006) the ques- 
tion then arises: Can radius be used to determine the mass of 
the star? 

Recent studies have shown interferometry has determined the 
radius of a star with a precision of about 1 % for the brightest 
stars and it has been suggested that we might obtain a precision 
of at least 4% for most near by solar-type stars whose diameters 
we could expect to nieasure jPiipers et alj2003l:lThevenin et al.1 
l2005t lDi Folco et al. 112004'). Our interest is exploring how this 
independent measure complements oscillation frequency infor- 
mation for solar-type stars. Can it determine the mass of the 
star? And if so, can we then use the frequencies to probe the 
physics of the star? How important is the radius for determin- 
ing other stellar parameters such as age and chemical composi- 
tion? What is the impact on the expected uncertainties? Now 
that we are expecting frequen cy errors less than about LO/j,Hz 
with the launch of CoRoT^ (Baglin et al.] l2000l) and Keplefl 
jBorucki et all 1 19971 1200 4: Basri"2004), does this have an ef- 
fect on our results? 

While the objective of this study is mainly to determine stel- 
lar mass, we also discuss results for other stellar parameters. 
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Section |2] describes the mathematical background, the physics 
of the models and the observables we use for this study. Sec- 
tion [3] discusses some theoretical results indicating the impor- 
tance of a radius measurement. Section |4] presents results of 
simulations which support the theoretical predictions. 

2. DESCRIPTION OF THE PROBLEM 

2.1. Mathematical Solution 



i 



The techn i que s incorporate d for this study follow that of 
ID (fl994) and ,Migho & MontalbanI (l2005l). For a 



rown et ; 

more e la borate descr i ption we refer the readers to lBrown et al.l 
(11994') or Press et al." ( 1992) and summarize here the main con- 
cepts for the purposes of understanding our work. 

2.1.1. Finding the parameters 

Taylor's Theorem allows us to approximate any differen- 
tiable function near a point Xq by a polynomial that depends 
only on the derivatives of the function at that point. To first 
order this can be written as 



f(x) = f(xo) + f'(xo)(x-xo) 



(1) 



where x = {-fj}^] are the parameters defining the system 
("input parameters") and f = are the M expected out- 

puts of the system that depend on x. These are the expected 
measurements or observables. To distinguish between the ex- 
pected observables and real observations we denote the latter 
byO = {0,}2i. 

We measure f and would like to find the set x that produce 
these measurements. This problem is a typical inverse problem 
whose solution x (in the pure linear case) falls neatly out from 

x = xo + VW-iu''5f, (2) 

where 

0-f(xo) 



Here 



5i-- 



f'(xo)e-i = D = UWV^ 



(3) 



(4) 



is the Singular Value Decomposition (SVD) of the derivative 
matrix (section |2.1.2| | and e = {e/jj^i are the measurement er- 
rors. 

In the case of a non-linear physical model, the method for 
finding the true values x is to use a modified version of equa- 
tion ^ iteratively. Incorporating a goodness-of-fit test such as 
a function 



2 ^ {m)-Oif 



(5) 



allows us to find the best set of parameters x by minimizing x^- 
There are many re liable algorithms available for performing 
such minimizations jPress et al.lll992h . The algorithm imple- 
mented in this work is the Levenberg-Marquardt algorithm be- 
cause it is known to be robust and incorporates derivative in- 
formation. It also usually converges within 2-3 iterations to its 
minimum. 

2.1.2. Singular Value Decomposition 

SVD is the factorization of any M xN matrix D into 3 com- 
ponents U, V'^ and W (eq. (|4]i). is the transpose of V which 
isanN xN orthogonal matrix that contains the input basis vec- 
tors for D, or the vectors associated with the parameter space. 



U is an M X N orthogonal matrix that contains the output ba- 
sis vectors for D, or the vectors associated with the observable 
space. W is a diagonal matrix that contains the singular values 
ofD. 

We use SVD in our analysis because it provides a method to 
investigate the information content of our observables and their 
impact on each of the parameters. Equations (O and (|4]l were 
defined in terms of the observed or expected error so SVD can 
be used to study various properties of D: 

• The expected uncertainties in each of the parameters via 
the covariance matrix 



(6) 



The significance of each of the observables for the de- 
termination of the parameter solution 



1/2 



(7) 



2.2. Parameters & Observables 



For this study we need to distinguish clearly between our pa- 
rameters and observables. The parameters are the input ingre- 
dients X to our physical model. For example in a stellar code the 
parameters include mass, age and chemical composition. The 
observables f are then the outputs of the models given these 
X and are those things that we can normally measure, such as 
magnitudes, effective temperature and oscillation frequencies. 
To refrain from ambiguity between the observable errors and 
the uncertainties in the parameters we shall denote the former 
by e and the latter by cr. 

2.2.1. Parameters & Models 

We describe a solar-type star with five adjustable parameters. 
These are mass M, age (or evolutionary stage) r, chemical com- 
position given by two of {X , Y, Z) where X + Y + Z= 1 and X, Y 
and Z are the initial hydrogen, helium and metal abundance re- 
spectively, and the mixing-length parameter a which describes 
convection in the outer envelope of the star. 

We use the Aarhus STellar Evolution Code (ASTEC) for 
stellar structure and evolution and the ADIabatic PuLSation 
code (ADIPLS) to calculate non-radial oscillation frequencies 
dChristensen-Dalsgaardl IT982I) . For each model evolved with 
the parameters {M,T,X,Z,a) ASTEC produces stellar quanti- 
ties which are then used to calculate the oscillation frequen- 
cies with ADIPLS. The outputs from the models implemented 
in this work are the global structure parameters such as radius 
effective temperature Teff, and a set of oscillation frequen- 
cies for / = 0, 1,2,3 and n = 1 ,2, ...,32. The Basel model atmo- 
spheres dLeieune et al.lll99"7h are used to calculate magnitudes 
and colors. R^„ T^ff, \ogg, and [Z/X] are used to obtain spectra 
which are then integrated in a range of wavelengths to produce 
magnitudes. 

We fixed the physics of the stellar evolution mod- 

els. The equation of state (EOS) is that of 

Eggle ton. Faulkner & Flannerv (1973). The o pacities are the 
OPAL 1995 tables (Jglesias & Rogers! Il996h . supplemented 
by Kurucz opacities at low temperatures. Co nvection is de- 
scribed by the classical mixing-length theory (iBohm-Vitensd 
1 195 8.) where the mixing-length £ is defined as £ = aHp, a is 
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the free mixing-length parameter and Hp is the local pres- 
sure scale height. Diffusion is not included and so we allow 
[M/H] ~ [Fe/H] to estimate {Z/X} where we are assuming a 
near solar-composition. We ignore overshoot effects and do not 
include Coulomb corrections. This study concentrates on one 
main-sequence model star with solar-type oscillations whose 
parameters are given in Table[T]but is extendable to a solar-type 
star with similar characteristics. 

2.2.2. Observables 

Typical observables for such a star come from various 
sources of observations. Spectroscopy provides effective tem- 
perature Teff, gravity \ogg and metallicity [M/H] measurements. 
With photometry we obtain magnitudes e.g. V and colors e.g. 
{U -V). Oscillation frequencies are obtained using time-series 
observations of photometric and/or spectroscopic origin. Inter- 
ferometry measures a limb-darkeneco stellar angular diameter 
in milliarcseconds (mas) and this coupled with a parallax gives 
an absolute diameter value (we use the radius Ri, value for this 
work). 

The observables (9, and expected errors e, for our model star 
are given in Table |2] The measurement errors reflect what is 
most currently quoted in the literature. However we do real- 
ize that some of these values are optimistic, and are unlikely to 
be improved upon. In this way, we can investigate how much 
further interferometry and asteroseismology can take us for de- 
termining stellar mass. 

We used frequency separations Afn.i and 6vaxi instead of in- 
dividual frequencies v^x where 

/S.Vn.\ = Vn.\-Vn-\.h (8) 
5l'n.0 = Vn+\X)-l'n.2- (9) 

While we realize that we are throwing away important infor- 
mation, the reasons for doing this are justified: the high sen- 
sitivity of the frequency values to changes in the parameters 
(non-linearity of the problem), and the unreliability of the treat- 
ment of the stellar surface layers and hence absolute values of 

We also assume that we can identify the (/,n) quantum num- 
bers. Rotation will induce a small effect (solar-type star) so we 
can assume (/,ot) ~ {l,m = 0). The degree / ca n be interpreted 
by us ing echelle-diagrams (see Figure 6 in iKieldsen et al.l 
(l2005h ). It is only the radial order n that is then difficult to iden- 
tify. This number is model-dependant, and there is no direct 
way of observing it. However with a rich frequency spectrum 
expected of these stars, the relative n position of each mode of 
degree / can be identified and then also the average frequency 
separations lS.v and 5v (independent of n). These observables 
are used for the initial fitting process, and an inspection of the 
i^n,! that come from the models allows an n identification to ±1. 

3. VARIOUS ROLES OF THE OBSERVABLES 

Why should a diameter measurement such as that expected 
from interferometry be important? This section highlights the 
importance that a radius measurement has when coupled with 
various observables and different measurement errors. 



3.1. Significance 

The observables play very different roles in the determination 
of the parameter solution depending on which combination of 
observables ((9,) are available and the size of their measurement 
errors (e,). Using equation (|7]i we calculated the significance 
of some observables 5(0,) using different observable combina- 
tions. 

Figure [T] (top panel) illustrates the importance of some typ- 
ical observables in the absence of a radius measurement and 
oscillation frequencies. The e, are those quoted from Table |2] 
While this figure is mainly for comparison with the lower panel 
we highlight a few points: the photometric observables appear 
to supply more information than the spectroscopic — this is 
contrary to what is often believed. For example, the colors 
are more significant than Teff. S([M/H]) ^ 0.85 where we may 
have assumed that most information about chemical composi- 
tion comes from [M/H] and so should be 1. Also note the 
high values of both S(logg) and SiV). They are responsible for 
determining both the radius and the mass of the star 

We then included R^,, two Ai^n j and one Sva.Q in the previ- 
ous set of observables. Figure [1] (lower panel) illustrates how 
the relative information content in each O,- changes. There are 
two notable changes: (1) \ogg and V contain almost no infor- 
mation while both the observed R^, and lS.Va\ become respon- 
sible for determining the true radius (and mass) of the star, (2) 
S(V -R) and S([M/H]) decrease by a small amount indicating 
that the information from colors and metallicity is not contained 
in the new observables. We believe this information is primar- 
ily the chemical composition. This implies that either colors or 
a metallicity are important observables to have when comple- 
menting seismic and radius measurements. We also note that 
while S(y-/?) decreases by > 25%, S>{U-V) remains the same. 

We included more than one Ai'n j to highlight how when 
we change the measurement errors, the information contained 
in one observable is passed to another, in particular, be- 
tween Ri, and Ai/n,i. If we had included just one Ai^n i, then 
S(Ai'n.i) ^ S((5i^n.()) '^l. When there are two observables con- 
tributing similar (but not the same) information, the content gets 
spread among them. Similarly, with three Ai/n,i, S(Ai'n,i) ~ 
0.55. Figure |2]illustrates this. S(Ai/n,i) (bold dashed line) and 
S(/?*) (bold continuous lines) are shown as a function of e{R) 
with eiv) = 0.5/iHz. At the smaflest e{R), S{R^)^\. Increasing 
e{R) causes S(/?*) to decrease, while S(At'n_i) increases because 
it is taking over its role. The non-bold lines show the same but 
for e{v) = 1.3,2.5/iHz. There is a similar but slower trend, be- 
cause the radius remains more important than Ai^n i even for 
higher values of e{R). Using only one Ai^n j would not show 
an increase in S(At'n,i) because it would be already near 1, but 
S(/?*) would decrease more slowly, while observables such as 
\ogg and V would increase in importance. 

Both Figures [T] and |2] show that the relative importance of 
observables depends not only on the individual e, but also on 
the combination of (9,- available, and that there is no straightfor- 
ward relationship between a set of O, and the parameters they 
are responsible for constraining. 



3.2. Parameter Correlations 



Information regarding the light intensity profile across the star is provided by the interferometric observations of the second lobe of the visibility function. In the 
case where the star is not fully resolved, i.e. we only get some visibility points on the first lobe, then we require a limb-darkened mod el for the star, rend ering 
model-dependent. The limb-darkened profile also varies if we use 1-D or 3-D atmospheric models. Both Aufdenb erg et al.l 120051) and lBigot et al.l j2006l) provide 
detailed and quantitative descriptions on this matter and also confirm the inadequacy of the 1-D models in reproducing these second-lobe measurements. 
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The combination of observables we use is important for the 
independence of the parameters. In general, there should be 
a unique set of frequencies for every possible combination 
of {M,T,X ,Z,a). However, we are not using individual i^n.i 
but frequency separations for the reasons explained in Sec- 
tion 12221 This creates some dependencies among the parame- 
ters such as that between M and t, M and a and M and X. With 
the addition of measurement errors, this degeneracy between 
different parameter combinations becomes worse. 

To illustrate some of the dependencies among the parame- 
ters, we calculated surfaces for two parameters p and q while 
keeping the other parameters fixed at their correct values. We 
evaluate 68.3%, 90% and 99% confidence levels as defined for 
a two-dimensional i. e. Ax^ = xlw, "Xp^ = [2-31, 4.61, 9.21]. 
We used e(R) = 0.01 R© and e(i^) = 1.3^Hz and defined x? = 

/i#observables — 2). 

Using a set of as observables the 99% confidence region 
encircled the correct parameter values to within less than 1 % of 
each of the corresponding parameters (it produced almost a dot 
centered on its true value), supporting the idea that each set of 
parameters should produce a unique set of ;/n,i. 

To show an example of the dependencies introduced by using 
the frequency separations, Figure[3]illustrates the correlated x^ 
surface for mass and a while holding (t, X, Z) fixed at their 
true values using OS3 (see section [33] ). Clearly for any value 
of M there is a corresponding a that will give a minimum x^ 
value. We hope that a minimization will lead to accurate val- 
ues, but we may need additional observables to help constrain 
the solution. 

3.3. Propagated errors 

Can we get a precise determination of the mass using a radius 
measurement? How well do we need to know our observables? 
Will the effort of obtaining an interferometric diameter be out- 
weighed by the benefits, or can some other observable produce 
a similar result? 

To answer these questions, we investigated the expected pa- 
rameter uncertainties through equation ^ for various observ- 
ables with different errors. For the remainder of this paper we 
shall discuss three sets of observables: 

OSl = {R,,T,f(,[M/HlAiy,Siy}, 
0S2 = {/;,,reff,[M/H],Ai/„,i} 

and 

OS3 = {R,, T,fi, [M/niAiy,^uSi^n.i} 

for n = 9, 10, . ..27 and / = 0, 1 , 2. For OS2 and OS3 we include a 
total of 21 large frequency separations, where we have assumed 
that we cannot identify every consecutive mode for each I and 
thus have even fewer frequency separations. For OS 3 we also 
include 5 small frequency separations. We remind readers that 
the parameter uncertainty is denoted by a while the measure- 
ment errors are denoted by e. 

For ease of reading we shall refer to e(i/) = 0.5/xHz as LN and 
1.3/iHz as HN. Figure |4] shows the theoretical uncertainties on 
each of the parameters as a function of e(R). We show results 
using LN (black) and HN (gray). The dashed lines show the re- 
sults for OS 1 . The continuous lines show OS2, and the dotted 
lines show OS 3. 

For all parameters we find that a{P) (parameter uncertainty) 
is a function of e(^) for all combinations of observables. We 



also see that for M, X and Z the different observable combi- 
nations and errors do not impact significantly the theoretical 
uncertainties for e(R) < 3%. Only for t and a, using 0S2 in- 
stead of OS 1 may not always be beneficial. This is because Sv 
is excluded from OS2 and this observable is essential for deter- 
mining r and, as it turns out, a too. 

• Mass: WTiatever combination we choose, we find that 
for e(R) < 3% a(M) does not change. A radius mea- 
surement will always be more important than seismic 
information once we reach a precision of this order As 
the radius error begins to increase, its relative weight for 
mass determination decreases, while the relative weight 
of the seismic information increases. The quantitative 
differences between using LN and HN for each observ- 
able combination for e{R) > 3% supports the idea that 
seismic information is the dominant contributor to the 
reduction of the uncertainty in mass at higher e(R). 

• Age: Either a reduction in frequency error or an increase 
in observables significantly affects o'(r). OSl and OS3 
do not benefit significantly from the radius measure- 
ment, but both produce better constraints on the age than 
OS2. This is precisely because they contain Siy which is 
fundamental in constraining this parameter However, 
for OS2 a precise radius measurement reduces (t{t) by 
200% for LN and HN. The lack of the observable Siy 
forces to constrain the age. Not surprisingly, we find 
that if we use LN instead of HN for all observable com- 
binations there is a large improvement in (t(t). 

• X: Adding observables and obtaining smaller frequency 
errors leads to subtle differences in <j(X) for low ra- 
dius error. OSl with either LN or HN produce very 
similar constraints on X. Using 0S2 does not change 
a{X) very much for e{R) < 3%. At larger values, we 
find that the additional large frequency separations (es- 
pecially for LN) constraint to about 10%, an improve- 
ment of a factor of 2 over HN. OS3H-HN and OS2H-LN 
produce similar results, indicating that it is not the pres- 
ence of Si' specifically that is important, just the addition 
of (seismic) observables. 

• Z: This is one of the parameters that is least influenced 
by either radius or seismic observables. It is mainly con- 
strained by [M/H], but because this implies knowing X, 
other observables play a role in constraining this param- 
eter Radius is most important for Z when we use OSl, 
because we have less information available. Using OS2 
and OS3 decreases the dependence of Z on In fact, 
with OS3H-LN, we find that the radius has almost no role 
in constraining Z. For Z, the radius is important only 
when there is not enough information available. As we 
increase the number of seismic observables, this infor- 
mation supercedes the role of the radius. We also find 
that whatever combination we use, a{Z) never decreases 
below 10% even using LN. 

• a: The determination of a seems to be more complex 
than the determination of any other parameter. In gen- 
eral, the addition of observables (in number) coupled 
with better frequency errors leads to lower parameter 
uncertainties. This is not always the case for a. OS2 is 
least influenced by a reduction in radius error In fact, 
we find for HN, OS 1 gives better constraints on a than 
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OS2 for e{R) < 4%. Only when an observable set con- 
tains the small frequency separation and for small val- 
ues of radius error, do we get a significant decrease in 
the theoretical uncertainty in a. Using LN instead of 
HN always affects the determination of a{a), but more 
so for OS2 when the radius information is less important 
than the frequency information. 

The role of the radius sometimes depends on what other in- 
formation is available. For example, if we seem to be missing 
information, then improving the radius error will lead to a cor- 
responding reduction in the uncertainty of some parameter, like 
Z. The determination of a shows an interesting trend. A precise 
radius determination complements precise frequencies. Either 
of these alone does improve the uncertainty, but both together 
have a larger impact on (T(a). Clearly for e{R) < 3% the choice 
of observable set does not change cr(M); it is which domi- 
nates the mass determination. 

3.4. Magnitude versus Radius 

Are these uncertainty values worth the observing time? We 
compared the propagated error using OS2 with My (magnitude) 
instead of R^, (radius) to see whether My could produce results 
similar to Ri,. My is the only other measurement that might 
constrain the mass as 7?^, does. 

Figure |5] shows a(M) as a function of both e{R) (continuous 
lines) and e{Mv) (dashed lines) for two e(7;ff) = [50,230/:] for 
LN (left panel) and HN (right panel). We have attempted to 
scale the expected errors on both observables so that the quoted 
errors are comparable in terms of measurement difficulty. We 
chose eiR) = 0.01 R© (- 1%) to be similar to e(Mv) = 0.05 
mag. For both cases cr(M) using My will never reach the pre- 
cision that a radius measurement allows. The smallest mass 
uncertainty using My is ^ 5%, while if we use Ri, it reaches be- 
low 1%. Also My needs to be complemented with other good 
observables such as T^ff and LN to get this precision, whereas 
Ri, is a completely independent measurement that does not need 
the help of other observables. These are two very clear reasons 
for choosing to obtain an interferometric measurement of di- 
ameter For HN (right panel) it is even clearer that 7?^, is most 
important. 

Even if our scaling is incorrect for the comparison of mag- 
nitude and radius errors, it is unlikely that an My measurement 
will have a precision of better than 0.05 mag, in which case the 
smallest mass uncertainty is really above 7%. 

4. SIMULATIONS 

To test the validity of the theoretical results presented in Sec- 
tion 13.31 we performed simulations and tried to fit the observ- 
ables to recover the input parameters. Part of this work also in- 
volved finding a reliable method to apply to any set of observ- 
ables so that we could succesfuUy estimate the parameters of 
the star within the theoretical uncertainties. We use two meth- 
ods: direct inversions and minimizations. The direct inversions 
provide a test for the linear approximation of equation ([T]l by 
using an initial guess of the parameters xq and equation We 
did this only for OSL Minimizations iteratively incorporate 
a modified version of equation |2](Levenberg-Marquardt algo- 
rithm) and were performed with OSl, OS2 and OS3. 

The observables O = {Rq, To, [M/H]o, f(n,i)o} were simulated 
as follows: 

Ro = Rr + r„e{R) (10) 



To = T, + r^e{T) (11) 
[M/H]„ = [M/H], + rge([M/H]) (12) 

Vio = Vk + rge(z^) + ( — ) (13) 

where R = R^,; T = T^ff; the subscipt 'o' denotes the simu- 
lated observed value; the subscript 'r' denotes the real value 
that comes directly from the stellar code; the f, are individ- 
ual frequencies; rg is a random Gaussian error for each ob- 
servable with (7=1; is an arbitrary reference frequency 
value, which was chosen to be 6000/iHz. Using these sim- 
ulated frequencies we calculated each of the Ai^n.i and Si^n.o 
from equations ([8]l and (|9]l. Table |2] provides the errors, and 
e(Ai'n.i) = e(Siy„,Q) = ^2e(t'n.i)^. 

4.1. Linear Approximation -vs- Minimizations for OSl 

Direct inversions will allow us to investigate the natural 
weight that Ri, and its error have on the determination of M, 
because we do not use a minimization algorithm but just the 
linear approximation and equation (|2]l. For each radius and fre- 
quency error, we generated a total of 100 simulations of the ob- 
servables. We obtained a list of xo from a grid of models (Sec- 
tion 14.2. It and chose the first set from this Ust to be the initial 
guess. For most of the simulations, an initial mass of ~ Q.9Mq 
was used (the real value = 1.03). The system is relatively linear 
in these parameters and observables hence the inversion results 
should be independent of the initial values. (This was tested, 
and the results varied to less than 1%.) The simulations were 
repeated for various frequency and radius errors. 

Figure |6] shows the mean and standard deviation of the in- 
verted {fit) mass from the simulations as a function of the ra- 
dius error. The dotted line shows the real mass value. Each 
panel corresponds to a different e{v). The upper panel shows 
LN, the center HN and the lower panel has eiv) = 2.5/iHz. The 
envelope dashed lines show the theoretical uncertainties. 

Figure |6] (top and center panels) shows that the mass uncer- 
tainty is a clear linear function of the error in radius. In fact, 
eiv) has a very modest effect on both the precision and the ac- 
curacy of the mass. For the two upper panels, there is no appar- 
ent difference, indicating that the radius primarily determines 
the mass. However at t{v) = 2.5/^Hz we do see an effect on 
both the precision ant/ the accuracy of the results. This appears 
to be a side effect of poorly constraining the other parameters, 
which indirectly leads to incorrect mass estimates. 

It is encouraging to find that the simulation results are con- 
sistent with theory and that we can safely estimate the mass to 
within its theoretical uncertainty. Unfortunately for e(R) = 5%, 
safely estimating the mass to within 20% may not be very use- 
ful. However, we do find that for €{R) < 2% we are capable 
of determining the mass better than the theoretical uncertainty 
and quite accurately (for the two smaller frequency errors). For 
all simulation results, the true mass value always falls within I 
standard deviation of the fit mass. 

For the other parameters, the theoretical predictions are con- 
sistent with the simulation results. Figure |7] shows the mean 
and the standard deviation of the inverted (fit) age. As expected, 
the size of the frequency error is very important, especially for 
OS 1 . e{R) does have a modest effect and is most noticable for 
LN (top panel). 

To test whether a minimization algorithm improves the un- 
certainties on the parameters, we repeated the simulations and 
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then fit for the parameters (using the minimizing method) for 
a range of radius errors and LN. Figure [8] shows the standard 
deviations of the fit parameters, (for age we show (t(t)/4). All 
of the parameter's real values were within the a shown in this 
figure. 

The precision of the parameters are better than the theoreti- 
cal ones (Figures l4l8] l. There is a clear dependence of the uncer- 
tainties on the radius error when t(K) < 3%. Unlike the theoret- 
ical uncertainties, for larger e{K) there is no significant degra- 
dation in the precision of the parameters. This uncertainty-error 
dependence changes drastically for larger t{K) because the min- 
imization method is forcing a balance between the errors and 
the observations, unlike when we use the simple linear approx- 
imation (Equation|2l). 

4.2. Simulations and Minimizations with 0S2 and 083 

Since we are introducing more frequency observables in OS2 
and OS3, we are effectively giving less weight to the radius 
measurement for constraining the parameters. The fit parame- 
ters are thus more sensitive to the initial parameter estimates. It 
is therefore necessary to use small grids of models to obtain a 
good initial guess, and then let the minimization algorithm fit 
aU of the observables. 

4.2.1. Model Grids 

We created grids of models that span the five dimensional 
model space over a range of values that encompass about 10- 
30% on either side of the parameter values given in Table [T] 
These values covered 0.95 < M < 1 .09 Mq, 0. 1 < t < 2.0 Gyr, 
0.68 <X < 0.76, 0.014 < Z < 0.022 and 1.2 < a < 1.8. We 
divided each of the parameter ranges into five (ten for mass), 
and calculated models for each of these points, and then in- 
terpolated only along the mass for a much finer grid. We are 
restricting the range in mass, age and Z, but we assume that we 
can constrain the parameters within this range from the observ- 
ables. 

We later extended the mass range to be between 0.75 and 
1 .20 Mq to investigate whether our original parameter range 
was too restricted. Using r = 0.1, X = 0.70, Z = 0.024, a = 1.5 
and varying mass, we used the grids to see how well we could 
recover the mass in this extended mass range. These parameters 
were chosen to test for any biases introduced by selecting our 
original ranges to be nearly symmetric about the correct values. 
We tested this for e{R) = 0.02 Rq and HN only. There were 100 
simulations for each mass. 

We searched through the grids to find all of the parameter 
combinations that reproduced and Aiy to within ±3e. We 
discarded the combinations with larger than 1, and adopted 
the mean value of each parameter from the remaining combi- 
nations to be the^f values. We later used them as the initial pa- 
rameter values for the minimization. Figure |9] shows the mean 
of the simulations ± 3(7 where a is the standard deviation of 
each group of simulations. The dotted line indicates the ideal 
position for the mean of the simulation results. While there 
seems to be some systematic offset as a function of mass, in 
most cases the 3<j crosses the correct value. For the other pa- 
rameters, the fit values are not as well determined (we haven't 
interpolated for them). This appears to be responsible for the 
offset in fit mass. 



4.2.2. Results of Minimizations 

" Extending the radius en'or to 0.001 Rq did produce a small improvement in the parameters at this level for LN, 



Initial estimates of the parameters were obtained from the 
model grids. We added a constant of 0.02 to the initial estimate 
of X for two reasons: l)X had always shown a systematic offset 
when using the model grids and, 2) if the value of from our 
initial estimates is too low, the minimization algorithm may not 
search among all possible solutions. We then allowed the code 
to minimize in all available observables. The solution usually 
converged within four iterations. We set the tolerance level 
to be 0. 1, where is calculated from equation (|5j, but we used 
3e errors for the process instead of le, and x^ = where n 
is the number of degrees of freedom. An ideal value of is 1 • 
Anything much larger implies a bad fit, an error in the physics, 
or incorrectly quoted errors. 

e(R) took values from O.OO1-O.O5R0 (- 0.1-6.0%). e{iy) 
took values of 0.5 (LN) and 1.3^Hz (HN). Figure [10] shows the 
mean and standard deviation (Icr) of the fit mass as a function 
of e{R) for 0S2H-HN. The dotted line is the true mass value. We 
show only this combination because those for and those 
using OS3 did not lead to a significant difference, except for t 
and a whose uncertainties are represented by the dashed lines 
in Figure nn 

Our global results show that for all of the parameters 1) the 
mean values fit to within Icr of the correct value (Figure [TOl 
shows this for M only). 2) The parameter uncertainties are 
smaller than those predicted by theory. This is highlighted in 
Figure [TT] 3) All of the uncertainties show some dependence 
on e(R). 

Mass is the parameter that is most affected by an improve- 
ment in the precision of Ri,. Figure [TO] shows the mean and 
standard deviation of the fit mass. a{M) is clearly dependent on 
e{R) for < 3% errors. If we reach 0.1% in radius error, then we 
can determine the mass to better than 1 %. For e(R) > 3%, a{M) 
does not change significantly. This can be seen more clearly 
in Figure [TT] For e(R)< 3%, there is little difference between 
the uncertainties arising from the minimizations (and more ob- 
servables) or the linear approximation with 5 observables (Fig- 
ures l6l8l ). This can be understood by considering that for small 
radius errors, this observable has such a large weight that it de- 
termines uniquely the mass value. In this regime, the different 
combinations of observables do not have any influence, just as 
the theoretical uncertainties showed (Figure [4] top panel). So 
whether we use a minimization method with 26 observables or 
a direct inversion with only 5 observables, having a well mea- 
sured radius will uniquely determine the mass. When the rela- 
tive weight of the radius decreases, i.e. its error increases, other 
observables begin to play a role. This is why we see a dramatic 
change in the mass uncertainty between using OSl, OS2 and 
OS3 when e(R) > 3% (Figures |ll0]l. 

Figure[TT]shows the standard deviations for all of the param- 
eters. The dashed lines represent the results for OS3H-LN and 
the dotted line for OS3H-HN (continuous are OS2H-HN). For M, 
X, and Z the results did not change significantly. 

• The top panel illustrates cr(T). The values are rather 
high, but the standard deviations that result from a star 
of 4 Gyr (instead of 1 Gyr) have similar absolute val- 
ues and thus much smaller relative values. Note the 
quaUtative differences between the set of observables 
with Si^ (dashed and dotted lines) and the set with- 
out (continuous line). Even if we use 21 large fre- 
quency separations with the precision of LN, we do 
not get the same age uncertainty as when we use OS3. 
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This demonstrates that 8v is crucial for obtaining a 
precise determinati on of the age (or st ellar evolution - 
ary stage) (e.g. Brown et al. (1994); Brown (|1994|) : 
iMigHo & Montalban (2005); Mazumdar (2005)). For 
e(^) < 2% we also see a slight drop in air). This is 
because for small radius errors, we get precise determi- 
nations of both M and X, so the range of possible ages 
that fit &v is restricted, and hence we get better age de- 
terminations. 

• The lower panel shows the standard deviations of all 
of the other parameters. Both Z and a exhibit a small 
e{R) dependence: 9-12% for Z and 7-10% for a over 
the radius error range. Theory (Figure HI predicted that 
a{Z) would never fall below 10%, and increase only 
slightly to 12-14% at larger radius errors. This trend 
is clearly consistent with the simulation results. (j{a) 
was predicted to be 20-22% over the range of radius er- 
ror for OS2H-HN and 8-14% for OS2H-LN. The simula- 
tions show no significant differences between these ob- 
servable combinations. For 0S3H-LN, theory predicted 
4-12% uncertainties, comparable to the dashed line in 
Figure [TT] a{a) benefits significantly from the addi- 
tional 5v reaching a precision of about 4% and slowly 
degrading with e{R) until e{R) = 3% when it levels off 
at a (T value comparable to the other observable combi- 
nations. (j{X) shows a dependence on radius error sim- 
ilar to (t(M). For e{R) < 3% we generally estimate X 
to within 2-3%. For larger radius errors, cr(X) does not 
increase, indicating the redundancy of a radius measure- 
ment after e{R) = 3%. 

While some parameters may not seem very dependent on 
e{R), having a radius observable allows us to estimate some pa- 
rameters correctly. This in turn leads to a better estimation of 
the other non- radius dependent parameters because of the more 
restricted acceptable range. 

These results stem from an automatic method of fitting the 
observations. Given that there is no human input, we are confi- 
dent in the results and expect that even more precise results will 
be possible by paying attention to each individual case, and in- 
dividual observable values. 

5. CONCLUSIONS 

We investigated the role that a precisely measured radius, 
such as that obtained through interferometry, has on the deter- 
mination of the mass of a star. We also looked at how radius and 
oscillation frequencies work together to determine other stellar 
parameters. 

• We found that the importance of a radius measurement 
depends on the combination of the available observables 
and their corresponding errors. For some typical ob- 
servables such as Ri, and Ai/„.i, we can expect a mass 
uncertainty of between 1 and 4% for e(R) between 0. 1 
and 3% (keeping in mind that these numbers stem from 
an automatic parameter search and in a real case study 
we would pay more attention to each observable and 



can expect better precision.) For e(R) < 3% the mass is 
uniquely determined by the radius observable, allowing 
us to use the frequencies just to probe the stellar interior. 

• Our simulation errors are consistent with the theoreti- 
cal errors when we use the linear approximation (Fig- 
ure |6]l. Using a minimization method yields smaller 
than predicted uncertainties with similar qualitative re- 
sults for all values of radius error, but particularly when 
e(^) > 3%. For these values, the weight of the radius 
measurement decreases, allowing other observables to 
play a greater role in the determination of the parame- 
ters. The minimization method allows an optimal trade- 
off between all observables and their corresponding er- 
rors. 

• We emphasize the importance of understanding how 
each observable contributes to the determination of each 
parameter. For example, if we want to get an estimate 
of Z better than 10% we need a more precise [M/H]. A 
radius measurement would not be useful. To correctly 
determine age (or stellar evolution stage) it is crucial to 
observe Sv. Even with 21 large frequency spacings, we 
will not determine the age to the precision that 1 small 
frequency spacing will provide. It is also interesting that 
Sv coupled with a small radius error had a strong impact 
on the precision of the mixing-length parameter a. 

• We know that a discrepancy exists between the observed 
and the model frequencies, which leads to a systematic 
offset in mass determination (e.g. Miglio & Montalban 
2005). By allowing the radius to determine the mass, 
we have the advantage that 1) we have an independent 
measurement that we can use to try to resolve the dis- 
crepancy, and 2) by using the radius to determine mass, 
we can use the frequencies to probe the stellar interior. 

• One final remark we would like to make is how effec- 
tive a radius measurement can be to detect some error in 
the model. During this study we conducted some hare 
and hound^tests. In one particular case, we found that 
all of the observations were fit very well except for the 
radius, which showed a larger than 3e deviation value. 
This turned out to be the result of an incorrectly quoted 
metallicity meausurement! This just shows how pow- 
erful the radius can be to detect an error/flaw, and this 
could possibly be used to detect a flaw in the physics of 
the models. 
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Table 1 



System Parameters 



Parameter 


Value 


M(Mq) 


1.030 


T (Gyr) 


1.00 


X 


0.740 


z 


0.018 


a 


1.50 



Table 2 
System Observables 



Measurement 


Value 


Error 




m 


(e,) 


(R©) 


0.946 




Teff (K) 


5421 


50 




4.5 


0.3 


[M/H] 


0.00 


0.05 


My (mag) 


4.50 


0.05 


(U-V) (mag) 


0.633 


0.005 


(V-R) (mag) 


0.400 


0.005 


Ai/ (^Hz) 


148.3 




5v (/iHz) 


14.6 




va,n (mHz) 


1996.9 





The errors on i?* and v are varied. 
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Fig. 3. — contour plot for M an d a. {X,Z,t') are fixed at their correct 
value. We used OS3 (see Section[33) with e(R)=0.02 Rq and e{u) = 1.3/iHz 
to calculate ■ 




5 10 15 ao 

e(R) (%) 



Fig. 2. — The change in the S(S*) and S(Ai'„ j) as we vary e(R). Af„ ; 
increases in significance as R^, becomes more poorly constrained. Bold lines 
indicate eiu) = 0.5/xHz, lighter fines indicate €(u) = 1.3, 2.5 fjHz. 
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Fig. 5. — Comparison of propagated mass error using either (continuous) or My (daslied) with 0S2 observables at two eCTeff) = [50,230ir|. The left panel is 
that when using ciu) = 0.5^Hz, the right is that for e(i/) = 1.3/uHz. The lower x-axis shows e(R) in appropriate units whereas the upper x-axis shows what we could 
consider a comparative t{Mv). 
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Fig. 6. — Mean and standard deviation of tlie fit (inverted) mass for the 
simulations as a function of e(S) using OS 1 . Upper Results for LN; Center 
Results for HN; Lower Results for e(i/) = 2.5;iHz. 



Fig. 4. — Comparison of the expected theoretical uncertainties between OSl 
(dashed), OS2 (continuous) and OS3 (dotted) observable sets and different 
frequency errors for all parameters as a function of e(R). Results for HN are 
shown by the gray hne, and those for LN are represented by the black lines 
(darker lines = smaller frequency error). 
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e{u) = O.b/xHz 



1.3/LiHz 



3 4 
e(R) (%) 




0.90 



0.80 



0.70 



0.70 0.80 0.90 1.00 

Real M (M^,) 



1.10 



1.20 



Fig. 9. — Mean and 3(t uncertainty of the fit mass as a function of the trae 
mass value resulting from the simulations using only the model grids. Dashed 
veilical lines indicate boundary of extrapolated and evaluated models. 



Fig. 7. — Similar to Figure|6] but for age. 
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Fig. 11. — Comparison of standard deviations of all of the fit parameters as 
a function of e(R). OS2+HN, OS3+LN and 0S3+HN (r only) are denoted by 
the soUd, dashed and dotted lines respectively. 



